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Abstract 

A recent study indicates that the a 2 a 2 s order QED processes of T — > J/ip + X decay are com- 
patible with those of QCD processes. However, in the endpoint region, the Non-relativistic QED 
(NRQED) calculation breaks down since the collinear degrees of freedom are missing under the 
framework of this effective theory. In this paper we apply the soft collinear effective theory (SCET) 
to study the color-singlet QED process at the kinematic limit. Within this approach we are able 
to sum the kinematic logarithms by running operators using the renormalization group equations 
of SCET, which will lead to a dramatic change in the momentum distribution near the endpoint 
and the spectrum shape consistent with the experimental results. 
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I. INTRODUCTION 



During the past 15 years, the interactions of non-relativistic heavy quarks inside quarko- 
nium have been understood to some extent using the framework of non-relativistic effective 
theories [l|, 0]. These theories reproduce the physics of full QCD or QED by adding local 
interactions that systematically incorporate relativistic corrections through any given order 
in the heavy quark velocity v. They provide generalized factorization theorems that include 
nonperturbative corrections to the color-singlet model, including color-octet decay mecha- 
nisms. All infrared divergences can be factored into nonperturbative matrix elements, so 
that infrared safe calculations of inclusive decay rates are possible 0] • These non-relativistic 
effective theories solve some important phenomenological problems in quarkonium physics. 
For instance, they provide the most convincing explanation to the surplus J/ip and ip' pro- 
duction at the Tevatron Q], in which a gluon fragments into a color-octet cc pair in a pointlike 
color-octet S-wave state which evolves nonperturbatively into the charmonium states plus 
light hadrons. The factorization formalism allows these fragmentation procedures to be 
factored into the product of short distance coefficients and long distance matrix elements 
among which the leading one is (O^^^Sx)) where are local four-fermon operators 

in terms of the non-relativistic fields. 

There are, however, some problems that remain to be solved. One challenging problem is 
with the polarization of J / ip at the Tevatron. The same mechanism that produces the J / ip 
described above predicts the J/ip should become transversely polarized as the transverse 
momentum p± becomes much larger than 2m c [5]. Though the theoretical prediction is 
consistent with the experimental data at intermediate p±, at the largest measured values of 
p± the J/ip is observed to be slightly longitudinally polarized and discrepancies at the 3cr 
level are seen in both prompt J/ip and ip' polarization measurements 0. 

A new problem arose as a result of measurements of the spectrum of J / ip produced in the 
T(IS') decay by the CLEO III detector at CESR [7j. NRQCD calculations have been made 
for the production of J / ip through both color-singlet and color-octet configurations [1, 0] . 
Theoretical calculations predict that the color-singlet process T(IS') — > J / ifjccg + X features 
a soft momentum spectrum. Meanwhile, the theoretical estimates based on color-octet 
contributions indicate that the momentum spectrum peaks near the kinematic endpoint 0] . 
In contrast to the theoretical predictoins, the experimentally measured momentum spectrum 
is significantly softer than predicted by the color-octet model and somewhat softer than the 
color-singlet case 0]. 

A more detailed study on the color-singlet contribution to this process has been presented 
recently [10!] . It was found that the contribution of the color-singlet QED process is compa- 
rable with the QCD process. NRQED calculations indicate that the QED process will give 
a large contribution to the spectrum near the end-point that is not observed in the data. 
This contribution results from the J/ip being produced back-to-back with a pair of gluons 
forming a low-mass jet. However, in this region of phase space, the NRQED calculation 
breaks down, since it does not contain the correct degrees of freedom. NRQED contains soft 
quarks, photons and gluons, but it does not contain quarks and gluons moving collinearly. 
The correct effective theory to use in situations where there is both soft and collinear physics 



is Soft-Collinear Effective Theory (SCET) [ll|, [l2j, [13|, \U 

A similar situation happens when studying e + e~ — > J/ip + X. The combination of 
SCET and NRQCD has been successful in reproducing the shape of the measured J/ip 
momentum spectrum in e + e~ — ► J/ip+X [15]. SCET has the power to describe the endpoint 
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regime by including the light energetic degrees of freedom. In addition, the renormalization 
group equations of SCET can be used to resum large perturbative logarithmatic correctoins. 
THE nonperturbative NRQCD martix elements arise naturely in deriving the factorizatoin 
theorem using SCET. 

In this paper, we use SCET to study the color-singlet contribution to the T — > J/ip + X 
decay near the endpoint via a virtual photon. We derive the factorization theorm in SCET 
for this process. We find that the spectrum is softer than the tree order prediction of 
NRQED when including perturbative and nonperturbative corrections near the endpoint, 
giving better agreement with the data than the previous predictions. 

II. MATCHING AND FACTORIZATION 

In this section, we derive the SCET factorization theorem for the color-singlet contribu- 
tion to T — > J/ip + X via a virtual photon near end-point. This factorization formula is 
crucial since the NRQED does not properly include the relevent collinear degrees of freedom 
and thus breaks down in this regime. This can be understood by analyzing the kinematics 
near the end-point. In the centre-of-mass (COM) frame, we have 



+ k x . (1) 

Here n = (1, 0, 0, 1) and n = (1, 0, 0,-1), we have defined z = (E^+p^) / My and r = m^/m^. 
We also assumed that = 2m c and My = 2rrib. k^ and k^ are the residual momentum of 
the QQ pair inside the T and J/ip respectively. Near the kinematic endpoint, the variable 
z — > 1 and thus the jet invariant mass approaches zero. In NRQED, an expansion of k^/mx 
is performed and hence the jet mode is integrated out, which is only valid when the jet 
mass is large compare to the residual momentum. The invariant mass of the jet is large 
away from the endpoint. As z — > 1, the jet becomes energetic, with small invariant mass. 
Hence we must keep k^/mx to all orders. As a result, the standard NRQED factorization 
breaks down at the endpoint. SCET is the appropriate framework for properly including 
the collinear modes needed in the endpoint in order to make reasonable predictions. 

To derive the factorization theorem in SCET, we start with the optical theorem in which 
the decay rate is written as 

2E ^=^^- b H I d'ye-^(T\O^y)\J/,p + X)(J/i, + X\O(0)\V, (2) 

where the summation includes integration over the X phase space, which includes both the 
ultrasoft (usoft) X u and collinear X c sectors. The SCET operator O is of the form 

O = J2 e-^+n>v C (jm, u>) V a ^ v J<#{ U ) O^O r , (3) 



Pi, 



Px 



2zMy + 2 H 
M- 



T 



(l--)n"+ (l-z)n» 
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where the Wilson coefficient C(/j,,u) is obtained by matching from QCD to SCET at some 
hard scale \i = The operator is contrained by the gauge invariance. In our case, to 
leading order we have 

J*(u) = Tr[5LXj> ( 4 ) 

& m = 4 (M-vfXc (5) 

<D» r = xl (A 2 -a)> b . (6) 

Here the A's boost the J/ip or T from the COM frame to an arbitrary frame, ip and x 
are the heavy quark and antiquark fields which create or annihilate the constituent heavy 
(anti-)quarks inside the quarkonia. The collinear gauge invariant field strength B" is built 
out of the collinear gauge field A" 

SJ = =iwJ {V a ± + g s (A a n , q )±)W n . (7) 
9s 



where 



ex P (~9spn ■ A n> A (8) 

perms ^ ' 



is the collinear Wilson line. The operator V is used to project out the large momentum 
label 13 . 

The hard coefficient containing the spin structure is obtained by matching the Feynman 
diagrams shown in fig. (U which gives 

_ , 327r 2 e c e b aa s ± ± 
La ^~ l 2N c m c m b 9a ? 9 ^ [J) 

where g^ u = g^ v — (n^n u + n u n^)/2. We have chosen the hard coefficient so that the Wilson 
coefficient C(fi,co) is 1 at the hard scale ji H . 

Inserting the operator in Eq. (jBJ) into Eq. (J2J), 0\y) picks an additional phase, and the 
differential rate becomes 

x {T\J«* 0^<D$(y)\J/^ + X) {J/4, + X\J^' & m O?(0)|T) 

= u)C{^ J) rV r aWv , A :tr aWv ' , (io) 
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In the exponent of Eq. (flQj) . we have used — Mrv^ + Vn tl /2 fa M T /2(1 — z)n M . Further- 
more, we can decouple the usoft modes from the collinear degrees of freedom using the field 
redefinition [14| 

K, q = YA^. (11) 

The fields with the superscript (0) do not interact with usoft degrees of freedom. In the 
color-singlet contribution the usoft Wilson lines Y cancel since Y^Y = 1. The T and the 
J/ip states contain no collinear quanta, so we can write 

A aPtxu,a'P'^'u' _ 1 / J4 -iM r / 2(1- z)n-y 

x (T\0^(y) a\a^ m Of (0)|T> (0\J^{y) J^"' (0)\0) . (12) 

Here we defined an interpolating field, a^, for the J/ip and used the completeness of states 
in the usoft and collinear fields 

^2\J/iP + X u )(J/i/; + X u \ = |J/^|=o^, (13) 

Y,\X C )(X C \ = I. (14) 

x c 

The T is a very compact bound state, due to the large b-quark mass. In a multipole 
expansion, long wavelength gluons interacts with the T color charge distribution through 
its color dipole moment since the state itself is color neutral. In the theoretical limit of very 



heavy bottom quark, this coupling to the dipole vanishes [181 ] . Therefore we are able to 
write 

AaPwafpfi'iJ _ 1 / j4 -iM r /2(l-z)n-y 

A "»' ~32^ b J dye 

x (T|G?(y) <#(0)|T> (0\O^(y) a^O^O) |0) (0\jf\y) J^' (0)\0) . (15) 

To proceed, we introduce the shape function for J/ip 

r d«- il+ (°l Xe?i ^c{y~) «W ^t a i Xc |0) 

s ^ )= JJL_ e -^- 1 __ J , (16) 

as well as the shape function for T 

Both shape functions are normalized so that J d/ + S^ t x(l + ) = 1. The are color-singlet shape 
functions can be related simply to the color-singlet NRQCD matrix elements 19| . 



( X ^^5(tn-d-k + )^a lX ) =e(k + )(^a^^a lX ), (18) 
which amounts to a shift from the partonic to hadronic endpoint. 
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In addition a jet function J L0 {k + ) is denned as 
(0|2> [BiBfi (y)Tr [B^] (0)|0) 



iV 2 



1 



a' 



dk+ 



The leading order result for the collinear jet function is 16l 

1 



in 



-0(/c H 



d^, 



£,(M r +Lu)/(2M r ) 



(19) 



(20) 



Using the spin symmetry relation 17 



A*Af(...a\..^...> 



-<^AfA? {...a k ... a k ...) 



(21) 



and applying the identity <? J Af A|' = 
J/V'j we can write the decay rate as 



(v 5 v 5 ' — g 55 '), where v s is the four- velocity of the T or 



dr 

dp^ 



r P[x,r) 



i 2 



|C(M T £,/x)| 2 e(M T - 2£ 



■X / 



in which 



9 N? 



1 + r 



(22) 



(23) 



and P[x,r] = (x 2 — 4r)(l + r)/(x(l — r) 2 ). Near the end-point, P[x, r] — > 1. The variable 
x is defined as x = E^jm^. It is straight forward to check that to the leading order the 
differential decay rate reproduces the NRQED calculation 10J . 



III. RESUMMING SUDAKOV LOGARITHMS AND PHENOMENOLOGY 



SCET has the power to sum logarithms using the renormalization group equations 
(RGEs). Large logarithms arise naturally in the processes involving several well-separated 
scales and will cause the perturbative expansion breaking down. By matching onto an ef- 
fective theory, the large scale is removed and replaced by a running scale \x. After matching 
at the high scale, the operators are run to the low scale using the RGEs. This sums all 
large logarithms into an overall factor, and any logarithms that arise in the perturbative 
expansion of the effective theory are of order one. 

In the previous section, we have matched onto the SCET color-singlet operator, by inter- 
grating out the large scale replacing it with a running scale \i. We now run the operator 
from the hard scale to the collinear scale, which sums all logarithms. The counterterm as 
well as the anomalous dimension used for running the operator in the RGEs have already 
been calculated in Ref. [16] , and we can lift the results from that paper. The result for the 
resummed differential decay rate is given by 



1 dr r 



P[x,r]Q(M r -2E X ) J dr] 



Ots{HH) 



i 27O7) 



(24) 
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where 7 is defined as 



2 

1 = Jo 



^(^ + h 2 + (l-^) 2 )( T ^ln, + iln(l- ?7 )))-^ 
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(25) 



To sum the large logarithms, the collinear scale \i c is chosen to be approximately m x and 



the hard scale is set to be \in = 2m c , in same way as in Ref. [lOf . 

The result from Eq. fl2~4"j) sums up the leading logarithmic corrections which are important 
only near the endpoint. Away from the endpoint, the logarithms that we have summed are 
not important and contributions that we neglected in the endpoint become dominant. We 
therefore would like to interpolate between the leading order NRQED color-singlet calcula- 
tion away from the endpoint and the resummed result in the endpoint. To do this, we define 
the interpolated differential rate as 

1 dr ( 1 drfA „. . \ 1 dr rcsum , , 

LO P[x,r] + (26) 



T dp^ \T dp^p J T dpip 

The first term in parentheses vanishes when approaching the kinematic limit, leaving only the 
resummed contribution in that region. Away from the endpoint the resummed contribution 
combines with the —P[x,r] to give higher order corrections in C( s (hh) to the spectrum. 

In fig. [21 we compare the resummed, interpolated decay rate, Eq. fl26|) . to the leading- 
order color-singlet result We use m c = 1.548GeV and m& = 4.73GeV. Aqcd is set to 
0.21GeV so that a s (2m c ) = 0.259. In our figure, the dashed line presents the leading-order 
color-singlet calculation and the solid curve corresponds to the interpolated decay rate with 
the collinear scale chosen as fi c = mx- The shaded band is obtained by varying the collinear 
scale from \i c = m^/y^ to jj c = \f2mx-, since the choice of scale could only be determined 
by higher order corrections. After resumming the spectrum shape softens near the end-point 
and is thus more consistent with experimental data Q. 



IV. CONCLUSION 

In this work, we study the color-singlet QED process for J / ip production in T decay in the 
kinematic limit region. Since the NRQED breaks down ate this limit, we apply the SCET 
to study the spectrum. Our calculation consists of matching onto a color-singlet operator 
in SCET by integrating out the hard scale. Once the usoft modes are decoupled from the 
collinear modes using a field redefinition, we are able to show a factorization theorem for 
the differential decay rate, in which the decay rate can be factorized into a hard piece, a 
collinear jet function, and usoft functions. As pointed out by Ref. 19| the usoft function 



in this case can be calculated, resulting in just a shift from the partonic to the physical 
endpoint. 

By running the resulting rate from the hard scale \xn to the collinear scale /i c , we sum 
the large Sudakov logarithms. Finally, we combine the SCET calculation with the leading 
order, color-singlet NRQED result to make a prediction for the color-singlet contribution 
via QED process to the differential decay rate spectrum over the entire allowed kinematic 
range. 
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FIG. 2: The decay rate 1/Todr/dp^ via QED process. The dashed curve is the tree level direct 
rate [10]. The solid line presents the interpolated resummed direct rate. The shaded band is 
obtained by varying the collinear scale from \i c = mx/V^ to /x c = \/2rnx, since the choice of scale 
could only be determined by a higher order calculation. 
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